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基于 改进 趾 传 算法 的 稀 踊 重 构 算 法 
潘 美 虹 +， 郑 ” JT 
(南京 航空 航天 大 学 电子 信息 工程 学 院 , 南京 211106) 


摘 XE: 和 实 疏 重 构 算 法 中 凸 松弛 法 在 恢复 效率 方面 、 贪 禁 追 踪 法 在 恢复 精度 方面 存在 不 足 ， 基 于 遗传 算法 和 迭代 优化 
的 思想 ， 结 合 模拟 退火 以 及 多 种 群 算法 的 优势 ， 提 出 了 基于 模拟 退火 遗传 算法 和 基于 多 种 群 遗 传 算 法 的 启发 式 稀疏 
重 构 算 法 。 所 提 算 法 均 从 传统 遗传 算法 易 陷 入 局 部 最 优 解 的 缺陷 出 发 ， 分 别 通过 保持 个 体 间 的 差异 性 和 提高 种 群 多 
样 性 来 搜索 待 求 稀疏 信号 的 全 局 最 优 解 ， 并 通过 理论 分 析 证 明了 所 提 算 法 参数 选取 及 搜索 策略 的 有 效 性 。 此 外 ， 以 
阵列 信号 处 理 中 空间 信 源 的 波 达 方向 (DOA) 估计 问题 为 例 ， 验 证 所 提 算 法 的 有 效 性 。 仿 真 结果 表明 ， 相 较 于 正 交 
匹配 追踪 OMP 算法 和 基于 11 范 数 奇异 值 分 解 的 11-SVD 算法 ， 所 提 算 法 提高 了 DOA 估计 的 精度 ， 且 降低 了 运算 复 
杂 度 ， 使 其 快速 收敛 至 全 局 最 优 解 。 
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Sparse reconstruction algorithm based on improved genetic algorithm 


q Pan Meihong", Zheng Qin 
> (College of Electronic & Information Engineering, Nanjing University of Aeronautics & Astronautics, Nanjing 211106, 
China) 


Abstract: Convex relaxation methods present drawback in terms of computational complexity, meanwhile, greedy pursuit 
methods have disadvantages in their reconstruction accuracy. Based on the inspiration of iterative optimization of genetic 
algorithm and combining with the advantages of simulated annealing and multi-population algorithm, we propose two 
heuristic sparse reconstruction algorithms based on simulated annealing genetic algorithm and multi-population genetic 
algorithm. The proposed algorithms aim at the defects of the traditional genetic algorithm that often trapped in the local 
optimal solutions. In order to solve this, we implement two strategies to search global optimal solutions of the sparse 
reconstruction via maintaining the differences among individuals and increasing the diversity of the population, respectively. 
The validity of the proposed algorithms on parameter selections and search strategy is proved by theoretical analysis. The 
proposed algorithms can be applied to the DOA estimation of multiple spatial sources in array signal processing to verify the 
effectiveness. Simulation results show that, compared with the OMP algorithm and 11-SVD algorithm, the proposed 


algorithms have improved the accuracy and reduced the computational complexity, which can converge to the global 


optimal solution in a fast manner. 


Key words: multi-population genetic algorithm; simulated annealing genetic algorithm; DOA estimation; sparse 


reconstruction 
0 引言 algorithm, GA) 是 一 种 求解 组 合 优化 问题 很 有 效 的 算法 。 目 
前 基于 遗传 算法 的 稀 政 重 构 方法 已 有 一 些 初步 的 研究 ， 文 献 


随 着 压缩 感知 理论 的 提出 ， 很 多 学 者 在 稀疏 重 构 方法 上 [4] 中 提出 了 一 种 模拟 退火 与 贪 禁 算 法 相 结合 的 稀疏 重 构 算 
做 了 大 量 研究 ， 提 出 了 很 多 有 效 的 重 构 算 法 ， 形 成 了 相对 完 法 ， 在 不 同 信 噪 比 下 对 图 像 进行 重 构 来 比较 该 算法 与 BP 算 
备 的 理论 和 算法 体系 ,并 在 处 理 实际 问题 中 得 到 了 广泛 应 用 。 法 以 及 OMP 算法 的 重 构 精 度 ， 结 果 表 明 相 同 条 件 下 ， 该 算 
目前 所 提出 的 重 构 算法 主要 分 为 贪 禁 追踪 法 和 凸 松弛 法 由 法 重 构 误差 更 低 。 文 献 [5] 提 出 了 一 种 基于 10 范 数 最 小 化 的 模 
两 类 。 这 两 类 算法 存在 各 自 的 优点 和 缺 隐 ， 贪 禁 算 法 在 恢复 ” 拟 退 火 稀疏 重 构 方法 ， 该 方法 相对 于 贪 禁 追踪 算法 ， 提 高 了 
效果 方面 不 如 凸 优 化 算法 ， 但 实现 起 来 相对 简单 快速 ;而 凸 重 构 精 度 以 及 重 构 效率 。 文 献 [6] 基 于 粒子 群 优 化 算法 ， 结 合 
优化 算法 的 计算 效率 不 如 贪 楚 算法 ， 但 具有 较 好 的 稳定 性 和 贪 禁 追踪 的 思想 ， 提 出 了 一 种 新 的 稀疏 重 构 算 法 ， 该 方法 比 
理论 保证 。 由 于 稀疏 重 构 问题 是 一 个 非 凸 的 组 合 优化 问题 ， 一 般 的 重 构 算法 具有 更 快 的 重 构 速 率 。 文献 [7] 把 粒子 群 优 化 
目前 还 没有 快速 精确 求解 该 问题 的 确定 性 方法 23。 而 前 面 提 算法 用 于 稀 玻 分 解 的 最 优 匹配 原子 的 搜索 来 降低 稀 玻 分解 复 
到 的 两 大 类 重 构 算 法 只 能 对 其 进行 近似 求解 ， 并 且 在 稀 玻 度 杂 度 ， 从 而 提高 用 稀疏 分 解 理论 进行 信号 处 理 的 计算 效率 。 
较 小 情况 下 不 能 获取 有 效 的 重 构 。 文献 [8] 将 遗传 算法 应 用 于 变换 域 稀疏 缺失 样本 信号 的 重 构 。 
稀疏 重 构 问 题 本 质 是 组 合 优化 问题 。 遗 传 算法 (genetic 文献 [9] 提 出 了 一 种 基于 遗传 算法 CGA) 的 测量 矩阵 优化 方 
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法 来 解决 压缩 感知 的 首 合 成 孔径 雷达 稀 玻 成 像 问题 。 文 献 [10] ” 计 ， 基 于 这 个 行 支撑 估计 ， 用 最 小 二 乘法 得 到 解 的 估计 : 
提出 了 一 种 基于 遗传 法 和 正则 化 先 验 模型 的 单 幅 图 像 SR -—— d 
重建 方法 ， 结 合 遗传 算法 和 迭代 收缩 算法 来 处 理 正则 化 先 验 ub d A 
模型 ， 避 免 陷 入 局 部 最 优 解 。 本 文 同时 采用 贪 禁 追踪 法 种 X X 是 该 稀 疏 重 构 的 问题 的 解 ， 则 可 以 得 到 
松弛 法 的 思想 ， 结 合 改进 的 遗传 算法 提出 了 新 的 稀 朴 重 构 算 | 44;Y -Y |l- min | AX -Yl T 
法 。 由 于 单一 的 遗传 算法 在 迭代 后 期 个 体 适应 度 趋 于 一 致 ， aud 
使 得 整个 种 群 的 进化 停滞 不 前 , 从 而 很 容易 陷入 局 部 最 优 解 。 上 述 式 子 等 价 于 
鉴于 模拟 退火 算法 更 强 的 全 局 搜索 能 力 ， 将 其 与 遗传 算法 相 — E: T 
结合 ， 在 降温 过 程 对 适应 度 进行 拉 伸 来 保证 算法 有 效 快速 地 
收敛 到 全 局 最 优 解 。 多 种 群 遗传 算法 则 突破 单 种 群 优化 的 框 式 (5) 和 (6) 说 明 稀 疏 重 构 的 本 质 就 是 寻找 符合 上 述 
架 ， 通 过 多 个 设 有 不 同 控制 参数 的 种 群 协同 进化 ， 同 时 兼顾 约束 的 行 支撑 估计 O 
了 算法 的 全 局 搜索 和 局 部 搜索 性 能 。 协 同 进化 过 程 中 保留 各 — “ ü i 
种 群 的 最 优 个 体 ， 并 通过 移民 算 子 将 最 优 个 体 传播 到 所 有 种 寻 此 可 以 定义 函数 : 
群 ， 和 迭代 结束 后 从 所 有 最 优 个 体 中 选 出 一 个 最 优 个 体 作为 最 ^, _ 1444Y -Y lle 
终 解 。 peg aet en 
本 文 将 所 提 算 法 用 于 阵列 信号 处 理 中 空间 信 源 的 DOA ”作为 遗传 算法 的 目标 函数 ， 即 适应 度 函 数 。 最 小 化 该 目标 函 
估计 来 验证 算法 的 有 效 性 。 目 前 基于 遗传 算法 的 DOA 估计 ” 数 来 得 到 稀疏 矩阵 X 的 行 支撑 估计 。 
方法 已 有 一 些 研究 : 文献 [11] 提 出 将 遗传 算法 应 用 于 MUSIC 优化 目标 : 
谱 峰 搜索 来 实现 DOA 估计 问题 ， 解 决 搜索 维度 过 高 而 导致 AO d 
的 算法 效率 低 的 问题 , 文献 [12] 则 结合 模拟 退火 与 遗传 算法 ， id 
将 其 应 用 于 加 权 子 空间 的 DOA 估计 中 来 提高 估计 精度 。 本 ”其 中 : @ 是 集合 3= 人 L2…, 刀 中 所 有 势 为 K 的 子 集 组 成 的 集 
文 将 稀疏 重 构 的 思想 与 改进 的 遗传 优化 算法 相 结 合 ， 用 于 合 。 
DOA 估计 , 并 与 正 交 匹配 追踪 算法 以 及 基于 奇异 值 分 解 的 11 2) 编码 方式 
范 数 最 小 化 的 凸 优化 算法 进行 比较 ， 验 证 算法 在 估计 精度 以 Xs Wi d m B dE o 行 所 在 位 置 集 
及 运算 效率 上 的 改善 。 ocsla-K.s-02..2), ERRER X 的 行 支撑 ， 作 为 种 群 
的 个 体 。 个 体 采 用 十 进 制 编码 ， 长 度 为 人， 个 体 上 的 基因 代 
1 ” 稀 玉 重 构 模 型 描述 表 非 零 行 所 在 的 位 置 标 号 。 
稀 玻 信号 x 经 过 感知 矩阵 A 压缩 测量 后 记 为 y=Ax ， 该 3) 初始 解 
方程 是 一 个 欠 定 方程 组 , 有 无 穷 多 个 解 , 但 由 于 x 的 稀疏 性 ， 设 种 群 规模 为 popsize， 每 个 解 表示 一 个 可 行 解 的 编码 。 
寻 此 通过 对 x 稀疏 度 的 最 小 化 约束 ， 方 程 就 可 以 得 到 唯一 的 。 ”为 了 让 算法 能 够 快速 收敛 到 全 局 最 优 解 ， 将 遗传 算法 初始 种 
解 。 考 虑 实际 中 更 一 般 的 情况 ， 建 立 含 噪 的 稀 疏 优化 模型 如 “和 群 中 的 个 体 设 定 为 比较 接近 最 优 解 的 支撑 集 集 合 ， 是 一 种 有 
T: 效 的 方法 。 本 文 遗传 算法 基于 贪 禁 追踪 的 思想 ， 采 用 阔 值 法 
t=arg min|xlh stlly-Axlb<e (1) ”来 设 定 初始 种 群 ， 得 到 一 个 初始 解 ， 然 后 在 该 初始 解 的 小 范 


H. rH 


tru xdé— MEL. x 中 大 多 数 元 素 为 0 或 者 接近 0; 


y 是 观测 向 量 y= 
相关 的 误差 界 ; 


元 素 个 数 ， 设 信号 稀疏 度 为 K， 则 天 才 xlh ， 满 
x 中 非 0 元 素 所 在 位 置 的 集合 


Ax+tmyeCm,xeCm，s 直 al 是 一 个 与 噪声 
灿 表 示 求 0 范 数 ， 用 来 求解 向 量 中 的 非 0 


HEK<M<«L, 


你 为 向 量 x 的 支撑 。 


考虑 更 一 般 的 情况 ， 观 测 信号 
Y-AX-«N,.YeC"T,AeC"-.XeC^T, NeQn, H 


采样 数 。 


] R(X)= 


示 和 矩阵 X 的 第 i 行 。 于 是 


下 : 


素 个 数 的 属性 ;| 


中 : 几 表 示 求 一 


号 为 时 域 的 多 次 采样 ， 记 为 


th: 了 表示 
üsisn|X,.2015& zs X 的 行文 撑 。 Xo 
建立 多 观测 信号 的 稀疏 重 构 模型 如 


JENE 


min |WX)|. st- AX -Y ||, có 
个 集合 的 势 ,集合 的 势 是 用 来 度量 身 
Ill- 表示 Frobenius 范 数 。 


st. 


(2) 
FH JU 


TL 
H 


2 ”基于 改进 遗传 算法 的 稀 琉 重 构 算法 


遗传 算法 包括 目标 函数 、 en 初始 解 、 
终止 条 件 五 个 基本 元 素 欠 
法 相 结合 ， 从 这 五 个 元 素 入 手 进行 模型 设计 。 


机 制 、 


1) 目标 函数 
Tr WEEK 


求解 该 问题 ， 


新 解 产 4 
将 遗传 算法 与 稀 琉 重 构 算 


FH 


己 知情 况 下 ， 问 题 式 (2) 可 以 转换 为 


min | AX -Y ||;  sz.|RQX)|e K 


首先 需要 寻找 一 


G) 
个 满足 | 人 上天 的 行 支撑 估 


围 


9, -(f | AY |; VP 


) 。 个 体 产生 方式 如 下 : 
(Q,; =Q; +r |r = randi(-10,10)) 


基于 以 上 理 


前 K 个 最 大 值 对 应 的 下 标 


集合 } 


内 进行 随机 波动 ， 来 得 到 遗传 算法 的 初始 解 集 合 
pop? = {0,0%,..., Oopsice 


(9) 


ib. Je H T BHUOB JOE SEDE DE Bs n d FJ 


算法 Csimulated annealing genetic algorithm based sparse 


reconstruction, SAGA-SR) 和 基于 多 种 群 
构 算 法 ( multi-population genetic algorithm based spars 
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pap 


reconstruction, MPGA-SR )。 从 单一 遗传 算法 易 


[避免 陷入 


解 的 缺陷 出 发 ， 通 过 改变 稀 玻 
于 全 局 最 优 解 的 概率 ， 从 而 减 小 算法 的 重 构 误 差 。 SAGA-SR 
利用 了 模拟 退火 算法 
局 部 最 优 解 ， 


以 一 定 概 率 接 受 较 差 解 的 
以 提 


MPGA-SR 则 


史 离 单 种群 优 化 的 框架 ， 


局 部 搜索 能 力 。 
和 专 算法 中 是 非常 重要 的 过 
更 新 


算法 的 计算 


量 和 搜索 结果 。 


不 同 算法 解 的 


此 本 文 针 对 两 种 算法 分 别 进行 


UA TD. 


2.1 


的 局 部 寻 优 能 力 ， 


SAGA-SR 算法 设计 
SAGA 结合 了 模拟 退火 算法 的 全 


使 得 求解 结果 收敛 于 全 局 


陷入 局 


XN 


CH (c 


部 最 人 4 


解 的 搜索 策略 来 提高 算法 收敛 


各 ， 从 


高 算法 的 全 局 搜索 能 力 。 


个 设 


制 参数 的 种 群 协同 进化 来 实现 最 优 解 的 搜索 。 该 算法 
硕 了 算法 的 全 局 搜索 和 所 
新 解 更 新 策略 在 遗 


ien 
机 制 不 


局 寻 优 能 力 和 遗传 算法 
最 优 解 的 


不 同 控 


同时 兼 


概率 提 
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gu 变异 概率 。 


1) 新 解 产 生机 制 
在 欠 代 优化 过 程 中 ， 为 了 保持 种 群 的 多 样 性 ， 引 入 一 个 
早熟 判断 标志 : 


At fi. c f (10) 
Ep: se 为 第 g 代 种 群 最 优 个 体 的 适应 度 值 ，; 
PEGATS; 了， 和 fr 分 别 表示 第 g VER, 
种 群 所 有 个 体 适应 度 均值 以 及 个 体 i 的 适应 度 值 ，H 表示 当 
前 种 群 中 个 体 适应 度 值 大 于 种 群 适应 度 均值 的 个 体 数 。 由 A 
动态 调节 交叉 和 变异 概率 : 


1 
P. = Pal 7 ome) (11) 
Egit 12 
Py = Puoll 1 ) (12) 


其 中 : pof pm 是 初始 的 交叉 和 变异 概率 。 
种 群 更 新 策略 基于 模拟 退火 的 思想 ， 产 生 下 一 代 种 群 : 


[e FSO) 
a4. r»exp-C(057)- f(05)/T) (3) 
其 中 ， i=l, popsie; 了 为 介 于 0 与 1 之 间 的 随机 数 。 
根据 种 群 中 个 体 的 值 ， 计 算 稀 朴 重 构 问 题 稀 疏解 : 
Ki =A. R: a -0 (14) 
2) 终止 条 件 
&E XLI EA SNR -1010g, AXE ， 给 定 一 个 最 优 解 的 


InA 
S, X JEVAROSNIRUMIELAB ER, BATF C140. 计算 得 到 ， 


对 该 最 优 解 S 求 目标 函数 的 上 界 。I 是 单位 矩阵 ，N 是 噪声 
矩阵。 


FOD IG 7 A&A;)N |l /AX +N |l; 
3I N le AY le 

3I N I. A AX" Il. - EN D) 

= qoem pa 


因此 可 以 令 e- 09 -D 作为 迭代 终止 的 条 件 。 从 算法 


(15) 


d) 计算 当前 种 群 个 体 的 适应 度 函 数值 f (Of) K =1,..., popsize , TR 
据 式 (13) 的 更 规则 ， 来 产生 下 一 代 种 群 。 

e) 更 新 温度 T=67 。 
终止 判断 : 当 (Ou) «e 或 8g>maxiter 算法 结束 ， 


得 到 最 优 解 


€, =arg min. Uf(QOY 。 将 该 个 体 代入 式 (14) 中 ， 即 可 输出 稀疏 


重 构 问题 的 最 终 解 外 i ; 否则 g=g+1， 转 a)~e)。 
2.2 MPGA-SR 算法 设计 
MPGA 利用 多 种 群 并 行 遗 传 思想 改进 遗传 算法 性 能 0 
不 同 种 群 赋 以 不 同 的 控制 参数 各 自 独立 进化 。 各 种 群 中 ， 交 
又 算 子 是 产生 新 个 体 的 主要 算 子 , 决定 算法 的 全 局 搜索 能 
变异 算 子 则 是 产生 新 个 体 的 辅助 算 子 ， 决 定 算法 的 局 部 搜索 
能 力 。 该 算法 能 够 避免 单 种 群 进化 过 程 中 出 现 的 过 早 收敛 现 
象 ， 并 且 也 能 加 快 收敛 速度 。 
1) 新 解 产生 机 制 
第 i 个 种 群 的 控制 参数 交叉 概率 和 变异 概率 ) 如 下 : 


i=1,..., popnum (16) 


ER: p.,p。e(0,D 是 常数 ， 分 别 表示 交 又 概率 的 下 限 以 及 变 
异 概 率 的 下 限 ; 7 和 4 分 别 是 介 于 0 与 1 之 间 的 随机 数 ; 
deam € (0.0.5) 分 别 表示 交叉 概率 以 及 变异 概率 最 大 可 变化 的 
幅度 ; popnum 表示 种 群 个 数 。 
选取 和 保留 每 个 子 种 群 的 优秀 个 体 ， 各 种 群 间 通过 移民 
算 子 进行 联系 来 促进 子 种 群 的 进化 ， 从 而 实现 多 种 群 协同 进 
化 。 移 民 操 作 如 下 : 
ar = Ooo 


Oor = arg max{ f (popi)} i=1,..., popnum (17) 
Ogs =arg min( f (pop:1)} 


第 i 个 种 群 中 个 体 的 更 新 策略 : 


rx I " f(Q < FOD 
Qf» fO) n 
Oe ASOD Uo 


b e 
x is KOLESO Dopnum 


的 计算 时 间 和 收敛 性 两 方面 考虑 ， 设 定 两 个 终止 条 件 ， 满 足 
任意 一 个 条 件 都 将 终止 算法 ， 即 : 如 果 种 群 中 最 优 个 体 的 适 


其 中 : Qee 为 第 i 个 种 群 第 a 和 个 个 体 进行 进化 操作 


应 度 函 数值 满足 fo.)<s， 或 者 迭代 次 数 达 到 参数 设 定 的 最 
大 允许 迭代 的 次 数 maxiter, 则 算法 终止 ; 当 条 件 S Oo) <e ^l 
满足 时 ， 可 以 提前 终止 算法 ， 避 人 免 多 余 的 迭代 ， 提 高 计算 效 
K, 当 算 法 陷入 局 部 最 优 解 时 ， 通 过 第 二 个 条 件 ， 可 在 有 限 
次 迭代 之 后 结束 ， 避 免 算法 陷入 死 循 环 。SAGA-SR 算法 的 
体 实现 步骤 如 下 所 示 。 

算法 1: 基于 模拟 退火 遗传 算法 的 稀疏 重 构 算 法 

输入 参数 : 感知 矩阵 A， 稀 玻 度 K， 种 群 规 模 popsize, 交叉 概率 Pr ， 
变异 概率  ， 初 始 温度 To, ， 冷 却 系数 纪 KIESA E, WKAR 
maxiter. 

初始 化 : Ek (1<K< popsize ) 个 个 体 编 码 为 Q , 
设 迭代 次 数 g-1, 4 T-Te. 

a) 选 择 。 依 据 目标 函数 f (OL). Kk =1,..., popsize 选择 最 优 个 体 保留 
下 来 ， 并 淘汰 最 差 的 个 体 。 
b) 交 叉 变异 。 以 概率 pc 选择 一 组 个 体 ， 个 体 间 随机 配对 ， 采 
单 点 交叉 的 方式 做 交叉 操作 ， 以 概率 pw 对 各 个 体 上 的 基因 做 变异 操 
作 ， 得 到 新 种 群 。 
c) 将 新 种 群 中 的 个 体 ,k=1,.…,popsize 代入 式 〈14) 得 到 
X, k =1,..., popsize ， 然 后 用 式 (7) 计算 当前 种 群 个 体 适应 度 fO) ， 
根据 式 〈16) 计算 早熟 判断 标志 As ， 然 后 由 式 (11) 和 (12) 更 新 交叉 


X (8) 确定 ， 


后 的 新 个 体 ， FO) 和 fO) 分 别 表示 遗传 操作 前 后 个 体 适 
应 度 值 ， 由 式 〈7) 计算 得 到 。 


算法 中 移民 算 子 的 功能 ， 以 源 种 群 中 的 最 优 个 体 代替 
标 种 群 中 的 最 差 个 体 ， 达 到 多 种 群 协同 进化 的 目的 。 同 时 选 
出 各 种 群 中 的 最 优 个 体 ， 并 将 其 放 入 精华 种 群 加 以 保存 ， 保 
证 各 种 群 产生 的 最 优 个 体 不 被 破坏 。 
2) 终止 条 件 
给 定 目 标 函 数值 的 上 界 < ， 计 算 方式 同 式 (15)， 从 精华 
种 群 中 选 出 最 优 个 体 Qw 。 若 该 个 体 满 足 f(Q%w)<s ， 则 算法 
终止 ， 当 算法 欠 代 次 数 达 到 最 大 友 代 次 数 仍 不 满足 目标 函数 
值 小 于 < 的 条 件 ， 则 算法 也 会 终止 。 设 置 最 大 迭代 次 数 是 为 
了 避免 求解 结果 陷入 局 部 最 优 解 时 ， 算 法 陷入 死 循 环 ， 作 用 
同上 。MPGA-SR 算法 的 具体 实现 步骤 如 下 所 示 。 
算法 2: 基于 多 种 群 遗 传 算法 的 稀疏 重 构 算法 
WASS: RAER A, AREK, FIR popsize， 子 种 群 数 
H popnum, FERRE i 的 交叉 概率 E, 变异 概率 pi, ,i=1， 
终止 参数 2 ， 最 大 迭代 次 数 maxiter。 
初始 化 : 产生 popnum 个 初始 种 群 pop4 ,k=1…, popnum ， 对 各 种 群 中 
个 体 进 行 编码 ， 编 码 方式 基于 式 (8), BARAKA g-1. 
a) 选 择 。 根 据 式 (7) 计算 所 有 个 体 目标 函数 值 
F(Q) i=1,..., popsize,Q, € pop, k =1,..., popnum ， 目 标 函 数值 越 小 个 体 


...popnum, 
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越 优 ， 从 每 个 子 种 群 中 选择 最 优 个 体 保留 下 来 ， 并 淘汰 最 差 的 个 体 。 
b) 对 popnum 个 子 种 群 进行 更 新 操作 : 
从 第 i 个 种 群 随机 选取 个 体 Q 和 Qj 按照 交叉 概率 pi 做 交叉 操 

作 产生 新 的 个 体 Ore W G, 计算 新 个 体 的 适应 度 值 , 利用 式 (18) 

对 每 个 种 群 中 的 个 体 进行 更 新 。 

对 交叉 操作 后 的 种 群 按照 变异 概率 pj, 

群 ， 处 理 方法 同上 。 

c) 对 popnum 个 子 种 群 进行 移民 操作 ， 设 置 初始 移民 次 数 m-0: 

基于 式 (17) 的 移民 策略 ， 找 出 第 站 个 种 群 的 最 优 个 体 Ou. DÀ 

及 目标 种 群 i+1 的 最 差 个 体 Qi ， 用 Qs BRE i+1 个 种 群 的 

Os 

判断 各 种 群 之 间 是 否 完成 移民 ， 若 m«popnum, W) m=m+1， 继 续 

移民 操作 ， 否 则 转 d). 
d) 从 每 个 子 种 群 中 选 出 各 自 的 最 优 个 体 ， 放 入 精华 种 群 中 保存 ， 

以 免 最 优 个 体 被 破坏 。 然 后 转 a)。 

终止 判断 : 在 精华 种 群 中 选 出 全 局 最 优 个 体 Qo ， 


进行 变异 操作 ， 并 更 新 种 
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为 了 7， 阵列 流 型 矩阵 为 4=[a(@),a(2)…,a(go)] ， 则 建立 信号 模 
型 : 


X(0) 2 AS()-N(D,t 91... T 
将 空间 角度 均匀 划分 为 L 等 份 ， 构 造 方向 矩阵 
A(9) 2[a4(0).a(0,)....a(0,) ， 其 中 : 


AN 


a(&) - [1 erdina. err aso y C FEL it EAS 7] 125 3E £1 R AE, 
建立 优化 模型 如 下 : 
X(t) 2 A(O)S, (0 - N(t,t =1,.. T 
in|X -AS,|, st|R(SOER 
对 算法 进行 有 效 性 验证 ， 以 估计 误差 作为 评判 准则 。 分 
别 对 以 下 两 种 情况 进行 蒙特 卡 洛 实验 , 每 一 点 进行 200 次 蒙 
特 卡 洛 实验 。 估 计 误 差 用 估计 值 与 实际 值 之 间 的 均 方 根 误 差 


来 描述 : 
1 
RMSE = EDS MEY (19) 


ES 
E 


OQ, arg, min (f(,,4)) ， 若 该 个 体 满 足 f(Q)<e 或 


g >maxiter ， 则 终止 算法 ， 得 到 最 优 解 ， 将 该 个 体 代 入 式 〈14) 中 ， 
即 可 输出 稀疏 重 构 问题 的 最 终 解 wm ;否则 g=g+1， 转 a)。 


3 ”算法 复杂 度 分 析 


SAGA-SR 算法 的 计算 量 主要 集中 在 初始 解 的 选择 过 程 
以 及 迭代 更 新 过 程 中 计算 目标 函数 的 过 程 ， 计 算 每 个 个 体 的 
初始 解 利用 了 冰 值 法 ， 该 步骤 的 计算 复杂 度 分 别 为 O(MLT). 
每 个 个 体 的 更 新 过 程 。 计 算 量 主要 集中 在 d) 计 算 目 标 函 数 的 
过 程 ， 其 复杂 度 为 OOQMKL2D)+OCMK7T)， 因 此 SAGA-SR 算法 
FRERE ZAE KX O(MLTxpopsize)+O(popsizex 
MK!®)+O(popsizexMKT) 。 考 虑 算法 的 最 大 迭代 次 数 为 
maxiter, 可 以 得 到 SAGA-SR 算法 复杂 度 为 O(MLTxpopsize)+ 
O(popsizexMK'! xmaxiter)4 O(popsizexMKT»xmaxiter) . 

MPGA-SR 算法 的 计算 量 主 要 集中 在 每 个 子 种 群 的 初始 
解 的 选择 、 和 迭代 更 新 过 程 目标 函数 的 计算 以 及 从 所 有 种 群 中 
选取 精华 种 群 这 三 个 过 程 中 。 考 虑 popnum 个 子 种 群 ， 每 个 
子 种 群 的 规模 为 popsize， 阔 值 法 选取 初始 解 ， 该 步骤 在 一 次 
迭代 中 的 计算 复杂 度 为 O(popnumxpopsizexMLT)， 算 法 aj 
算 个 体 bs B XQ € 48 HE 为 O(popnumx 
popsizexMKÜ?!)« O(popnumxpopsizexMKT), d)3 BUR *E qb E 
计算 复杂 度 为 O(popnumxMK 0?) OpopnumxMKT)。 针 对 整 
个 算法 的 最 大 迭代 次 数 maxiter, 可 以 得 到 MPGA-SR 算法 复 
林 O(MLTxpopsizexpopnum) 


[2 


4O(popsizex — (popnum-1)x 


MK"! xmaxiter)4O(popsizex (popnum+1) x MKTxmaxiter) 。 
4 ”算法 的 有 效 性 验证 

贪 森 追踪 类 算法 选用 典型 的 OMP 算法 [9， 西 优化 算法 
采用 奇异 值 分 解 h 范 数 最 小 化 算法 器 ， 对 OMP, n-SVD. 
GA-SR、SAGA-SR、MPGA-SR 这 五 种 算法 进行 性 能 比较 。 
4.1 DOA 估计 问题 

本 节 针 对 天 - 稀 疏 信号 的 恢复 问题 , 以 空间 信 源 的 波 达 方 
向 (direction of arrival) 估 计 09 作 为 仿真 实验 ， 将 本 文 提 出 算 
法 与 贪 禁 追 踪 类 算法 中 的 正 交 匹配 追踪 算法 (OMP)、 凸 松 
弛 类 算法 中 的 基于 1 范 数 约束 的 凸 优化 算法 (11-SVD) 法 进 
行 比 较 ， 以 验证 本 文 所 提 方 法 在 估计 精度 和 估计 效率 上 的 性 
能 提高 。 

考虑 K 个 窄带 相干 信号 分 别 从 9,2…,Gk 入 射 到 M 个 各 
向 同性 的 传感器 构成 的 均匀 线 阵 上 。 阵 元 间距 为 d， 快 拍 数 


其 中 : O 为 蒙特 卡 洛 实验 次 数 ，K 为 信 源 个 数 ，&6, 表示 第 gq 


次 蒙特 卡 洛 实验 对 第 个 信号 的 估计 值 ; 
0, =22.6°,0, =27.2°,0,=35.9° 表示 第 大 个 信号 的 真实 值 。 
4.2. 重 构 误 差 随 快 拍 数 及 信 噪 比 的 变化 

考虑 一 个 阵 元 数 为 m=20 的 均匀 直线 阵列 。 智 能 算法 参 
数 设置 : 种 群 规 模 为 100, 交叉 概率 为 0.85, 变异 概率 为 0.05， 
进化 代数 为 50。 模 拟 退 火 遗 传 算法 的 初始 温度 为 5， 降 温 系 
数 为 0.9。 多 种 群 遗 传 算法 ， 种 群 数 为 5， 种群 规 模 为 100， 
各 种 群 的 交叉 概率 是 0.8-0.9 间 的 随机 数 ， 变 异 概率 为 
0.05~0.1 之 间 的 随机 数 。 

仿真 1 假设 入 射 信号 为 3 个 远 场 窄带 独立 信号 , 信 源 的 
AG fü BE Ay 3) Jg 6 =22.6°,0, =27.2°,0, =35.9° 。 固 定 信 噪 比 为 
10 dB， 快 拍 数 从 100 开始 以 20 为 步 长 变化 到 400, K 1(a) 
为 五 种 算法 估计 结果 的 均 方 根 误 差 随 快 拍 数 变化 曲线 。 

仿真 2 假设 入 射 信号 为 3 个 远 场 窄带 独立 信号 , 信 源 的 
入 射 角度 分 别 为 和 =22.6".2 =27.2°,0, =35.9 。 固 定 快 拍 数 为 
200， 信 噪 比 从 -10 dB 开始 以 2 dB 为 步 长 变化 到 2 OdB, 
1(b) 为 五 种 算法 估计 结果 的 均 方 根 误差 随 信 品 比 变化 曲线 。 

1 在 不 同 快 拍 以 及 不 同 信 噪 比 下 的 仿真 结果 表明 所 提 
算法 的 DOA 估计 误差 明显 优 于 OMP 算法 。 图 1(a) 的 实验 结 
果 表 明 MPGA-SR 算法 比 h-SVD 算法 估计 性 能 更 优 。 图 1(b) 
表明 在 信 噪 比 低 情 况 下 ，1-SVD 算法 估计 误差 明显 大 于 其 他 
算法 ， 表 明了 本 文 所 提 的 两 种 算法 在 低 信 噪 比 下 估计 的 有 效 
性 。 对 比 三 种 基于 遗传 算法 的 估计 结果 ，MPGA-SR Sf 
估计 精度 优 于 SAGA-SR, H. SAGA-SR 算法 的 估计 精度 优 于 
GA-SR 算法 ， 验 证 了 改进 遗传 算法 的 有 效 性 。 

4.3” 重 构 误 差 随 角 度 间隔 及 稀 下 度 的 变化 

算法 参数 与 5.2 节 设 置 相同 。 阵 元 数 设 为 20， 固 定 信 噪 
HE 10 dB, 快 拍 数 200。 对 算法 的 角度 分 辩 性 能 以 及 不 同 稀 玻 
度 下 的 估计 误差 进行 仿真 。 

仿真 3 固定 其 中 一 个 入 射 角度 为 35.7*?， 另 外 一 个 角度 
与 固定 角 的 之 间 的 间隔 从 1° 开 始 ， 以 1° 为 步 长 变化 到 309. 
图 2(a) 为 五 种 算法 角度 偏差 随 空间 入 射 角 之 间 的 角度 间隔 变 
化 曲线 。 

仿真 4 空间 信 源 的 入 射 角度 为 9=22.6* ， 多 =27.2? ， 
4 =35.9" 。 通 过 变化 阵 元 数 来 改变 观测 信号 的 稀疏 度 。 稀 玻 
度 定义 为 KM，K 为 信 源 数 ，M 为 阵 元 数 。 阵 元 数 从 3 开始 
以 3 为 步 长 变化 到 36. 图 2(b) 为 五 种 算法 估计 结果 的 均 方 根 
误差 随 稀疏 度 变化 曲线 。 
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显 大 于 其 他 


I 种 算法 ; 
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